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A  theoretical  modeling  approach  is  presented,  which  describes  the  behavior  of  a  typical  fuel  cell-heat 
engine  hybrid  system  in  steady-state  operating  condition  based  on  an  existing  solid  oxide  fuel  cell  model, 
to  provide  useful  fundamental  design  characteristics  as  well  as  potential  critical  problems.  The  different 
sources  of  irreversible  losses,  such  as  the  electrochemical  reaction,  electric  resistances,  finite-rate  heat 
transfer  between  the  fuel  cell  and  the  heat  engine,  and  heat-leak  from  the  fuel  cell  to  the  environment  are 
specified  and  investigated.  Energy  and  entropy  analyses  are  used  to  indicate  the  multi-irreversible  losses 
and  to  assess  the  work  potentials  of  the  hybrid  system.  Expressions  for  the  power  output  and  efficiency 
of  the  hybrid  system  are  derived  and  the  performance  characteristics  of  the  system  are  presented  and 
discussed  in  detail.  The  effects  of  the  design  parameters  and  operating  conditions  on  the  system  perfor¬ 
mance  are  studied  numerically.  It  is  found  that  there  exist  certain  optimum  criteria  for  some  important 
parameters.  The  results  obtained  here  may  provide  a  theoretical  basis  for  both  the  optimal  design  and 
operation  of  real  fuel  cell-heat  engine  hybrid  systems.  This  new  approach  can  be  easily  extended  to  other 
fuel  cell  hybrid  systems  to  develop  irreversible  models  suitable  for  the  investigation  and  optimization  of 
similar  energy  conversion  settings  and  electrochemistry  systems. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Nowadays  the  development  of  clean  energy  systems,  for 
both  transportation  and  stationary  applications,  is  recognized  as 
mandatory  to  satisfy  well-known  environmental  and  regulatory 
requirements  in  terms  of  emissions  and  energy  conversion  effi¬ 
ciencies.  Because  of  their  high  efficiencies  and  zero  toxic  emission 
levels,  fuel  cell  systems  are  considered  as  one  of  the  most  attrac¬ 
tive  solutions  in  the  automotive  and  power  generation  industry. 
Among  the  existing  fuel  cell  technologies,  the  high-temperature 
fuel  cells  show  a  great  promise  due  to  their  high  energy  con¬ 
version  efficiency,  fuel  flexibility  and  high  temperature  of  the 
exhaust  heat  which  allows  favorable  co-generation  and  combina¬ 
tion  with  other  types  of  power  generators  such  as  heat  engines 
[1-4]. 

After  the  concept  of  fuel  cell-heat  engine  hybrid  systems  was 
proposed  at  early  1990,  a  number  of  investigations  have  been 
carried  out  about  the  system  performance  from  the  energetic 
point  of  view  [5-7].  In  a  group  of  works  [8-10],  the  cell  perfor¬ 
mance  was  evaluated  based  on  the  operating  curve  interpolated 
from  experimental  test  results.  Some  authors  [5,6,11]  performed 
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a  parametric  analysis  to  understand  the  effects  of  different  oper¬ 
ating  conditions  on  the  performance  of  a  specified  system.  Most 
of  the  papers  are  dealing  with  the  theoretical  cycle  analysis  and 
simulation  of  the  possible  configurations  of  the  hybrid  systems 
[11,12]. 

Since  it  is  difficult  to  experimentally  quantify  the  interrelated 
parameters  governing  a  hybrid  system,  theoretical  modeling  and 
numerical  analysis  become  essential  for  the  optimization  of  the 
system  design  and  operating  conditions.  Therefore,  the  purpose  of 
the  present  paper  is  to  clarify  these  issues  and  present  a  general 
and  fundamental  analysis  of  the  theoretical  performance  potential 
of  a  fuel  cell-heat  engine  hybrid  system. 

Based  on  the  simplifying  assumptions  derived  from  literature 
and  performance  data  provided  in  [13],  an  irreversible  thermody¬ 
namic  modeling  for  a  fuel  cell-heat  engine  hybrid  system  is  further 
carried  out  to  simulate  the  system  response  to  the  changes  in  the 
operating  conditions  and  design  parameters  through  parametric 
study.  The  different  sources  of  energy  losses  are  specified  and  the 
performance  characteristics  of  the  system  are  investigated.  In  par¬ 
ticular,  both  the  heat-leak  from  the  fuel  cell  to  the  environment 
and  heat  transfer  between  the  fuel  cell  and  the  heat  engine  are 
considered.  As  a  consequence,  the  performance  of  the  hybrid  sys¬ 
tem  is  optimized.  These  findings  indicate  which  parts  of  the  system 
deserve  the  greatest  effort  if  one  wants  to  improve  the  overall  cycle 
efficiency. 
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2.  An  irreversible  model  of  the  fuel  cell-heat  engine  hybrid 
systems 

Although  fuel  cell  technology  has  been  studied  extensively,  the 
best  way  to  employ  fuel  cell  unites  for  the  generation  of  electrical 
power  remains  to  be  determined.  Optimal  thermal  management 
allows  for  effective  use  of  the  system’s  byproduct,  heat,  leading 
to  substantial  increases  in  the  overall  system  efficiency.  The  heat 
generated  in  a  fuel  cell  may  be  dissipated  by  convection,  conduc¬ 
tion  or  radiation.  However,  mostly  the  heat  may  be  used  to  run  a 
thermodynamic  cycle,  such  as  a  heat  engine,  for  additional  power 
generation. 

It  is  our  opinion  that  initially  a  fuel  cell-heat  engine  hybrid  sys¬ 
tem  has  to  be  as  simple  as  possible  to  achieve  competitive  first 
cost,  and,  probably  even  more  important,  high  reliability  and  avail¬ 
ability.  Thus,  an  illustration  of  a  simplified  fuel  cell-heat  engine 
hybrid  system  is  shown  in  Fig.  1 ,  where  fuel  cell  acts  as  the  high- 
temperature  heat  reservoir  of  a  heat  engine  for  a  further  production 
of  power.  The  role  of  the  regenerator  in  the  hybrid  system  is 
to  preheat  the  incoming  fuel  and  air  with  the  high-temperature 
exhaust  gas  of  the  fuel  cell  and  to  ensure  that  the  fuel  cell  works 
at  steady-state.  This  model  is  helpful  to  understand  the  architec¬ 
ture  of  any  combined  cycle  of  fuel  cells  and  heat  engines  and 
can  be  easily  expanded  to  analyze  different  chemical  reaction 
processes. 

The  whole  hybrid  system  is  formulated  based  on  the  following 
general  simplifications  and  assumptions  [13-18]: 


1.  Both  the  fuel  cell  and  heat  engine  are  assumed  to  operate  under 
steady-state  conditions. 

2.  All  gases  are  assumed  to  be  ideal. 

3.  Perfect  regeneration  is  considered  for  the  heat  exchanger. 

4.  Operating  temperature  and  pressure  are  uniform  and  constant 
for  the  fuel  cell  domain  including  the  inlet  reactants  and  outlet 
products. 

5.  Complete  chemical  reactions  are  considered  and  no  reactants  are 
remained  after  the  reaction. 


With  the  help  of  these  assumptions,  the  governing  equations  of 
this  hybrid  model  will  be  presented  in  this  section  as  the  influence 
of  some  new  parameters  is  considered. 


▲ 

I  Pcell 


Input  Output 


Fig.  1.  The  schematic  diagram  of  a  typical  fuel  cell-heat  engine  hybrid  system. 


2.1.  A  fuel  cell  model 

The  fuel  cell  model  presented  in  this  study  has  been  reported 
in  Ref.  [13],  which  is  based  on  a  solid  oxide  fuel  cell  (SOFC)  using 
hydrogen  as  fuel  and  air  as  oxidant.  The  basic  working  principle  is 
illustrated  in  Fig.  2.  The  overall  electrochemical  reaction  in  such  a 
cell  can  be  summarized  as  H2  +  (l/2)02  -*  H20  + Heat  +  Electricity. 

It  is  well-known  that  the  measured  open-circuit  potential  in  a 
practical  fuel  cell  is  usually  lower  than  the  ideal  reversible  poten¬ 
tial,  because  there  always  exist  some  irreversible  losses  originating 
primarily  from  activation  overpotential  (Vact),  ohmic  overpoten¬ 
tial  (V0hmX  and  concentration  overpotential  (VCOnc)  [16,19-24]. 
One  way  to  handle  this  behavior  is  to  assume  some  electronic 
current  leakage  through  the  electrolyte  [21,25-29].  Thus,  a  leak¬ 
age  resistance  can  be  introduced,  which  is  in  parallel  with  the 
external  load.  When  electrons  flow  through  the  internal  and  exter¬ 
nal  fuel  cell  circuit,  the  local  behavior  of  the  irreversible  losses 
is  described  as  a  simple  equivalent  electrical  circuit,  which  is 
series  of  a  reversible  voltage  determined  by  the  Nernst  equation, 
an  internal  resistance  made  up  by  the  sum  of  three  overpoten¬ 
tial  contributions,  and  a  leakage  resistance  in  parallel  with  the 
load. 

The  maximum  electrical  work  obtainable  in  a  fuel  cell,  operat¬ 
ing  at  constant  temperature  (T)  and  pressure  (p0  =  1  atm),  is  given 
by  the  change  in  Gibb’s  free  energy  (-AG)  of  the  electrochemical 
reaction.  The  relationship  between  the  Gibb’s  free  energy  (G)  and 
the  reaction  enthalpy  (H)  is  known  to  be 

-AH  =  -AG  -T  AS  (1) 

where  Q.  =  -  AH  corresponds  to  the  total  thermal  energy  avail¬ 
able  in  the  system,  AS  denotes  the  change  in  entropy  and  -TAS 
represents  the  amount  of  heat  generated  by  a  fuel  cell  operat¬ 
ing  reversibly.  The  enthalpy  change  and  Gibbs  free  energy  change 
between  the  products  and  the  reactants  of  the  global  electrochem¬ 
ical  reaction  at  temperature  T  for  the  steady-state  fuel  cell  can  be, 
respectively,  expressed  as 

AH=-^Ah  (2) 


and 


AC=-^Ag(T,p) 

Her 


(3) 
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where  Ah  =  -  Yvjhj,  A g(T,  p)  =  ^vkpik(T,  p)  - 

k  j  k 

JfvjUjiCp),  ne  is  the  number  of  electrons  transferred  in  reac- 

j 

tion,  i  is  the  current  density,  A  is  the  surface  area  of  the  polar 
plate  (supposing  that  the  bipolar  plates  have  the  same  area), 
F=  96,485  C  mol-1  is  Faraday’s  constant,  h  is  the  molar  enthalpy 
of  the  species,  fi  is  the  partial  molar  Gibbs  free  energy  of  species 
(which  is  also  the  chemical  potential),  v  is  the  stoichiometric 
coefficient  for  species,  and  subscripts  k  and  j  represent  the  fcth 
product  and  jth  reactant  of  the  reaction,  respectively.  Specif- 
ically,  Ag(T,p)  =  Ag°(r)-Rrin(pH2Po/22/PH2o)  denotes  the 
molar  Gibbs  free  energy  change  for  the  fuel  cell  reaction,  where 
R  =  8.314  J  mol-1  K-1  is  the  universal  gas  constant,  pH2,  Po2»  and 
Ph2o  are  the  partial  pressures  of  reactants  H2, 02,  and  H20,  respec¬ 
tively.  It  is  noteworthy  that  Ag°{T)  =  Ah° -T  As°  is  the  molar 
Gibbs  free  energy  change  at  p0  =  1  atm  and  is,  therefore,  called  the 
standard  molar  Gibbs  free  energy  change,  which  also  depends  on 
temperature  [14,16,19,23,30,31  ],  and  that  the  calculation  of  Ag°(T) 
is  based  on  the  tabulated  values  [32]  at  operating  temperature  T. 

Since  the  Gibbs  free  energy  change  of  an  electrochemical  reac¬ 
tion  is  a  measure  of  the  maximum  electrical  energy  obtainable 
from  the  reaction  [30],  the  maximum  (reversible)  power  output 
generated  by  the  reaction  in  the  fuel  cell  can  be  determined  by 


j  a 

Prev  =  -AG  =  -—=Ag(T,  p)  (4) 

Ii01 

As  mentioned  previously,  Eq.  (4)  is  a  measure  of  the  maximum 
work  obtainable  from  the  reaction  in  the  reversible  fuel  cell.  In 
practice,  this  part  of  energy  is  never  completely  utilized  because 
of  the  various  thermodynamic  and  electrochemical  irreversibili¬ 
ties  [33,34].  When  the  fuel  cell  works  normally  and  produces  useful 
power  through  the  external  load,  combining  those  irreversibilities 
yields  the  rate  of  the  total  entropy  production  of  an  irreversible  fuel 
cell,  which  includes  the  entropy  production  rate  resulting  from  the 
internal  resistance  and  leakage  resistance  as 


$tot 


^int^nt  ^leak^ea^ 

To  To 


(5) 


where  T0  is  the  ambient  temperature,  Iint  and  Iieak  represent  the 
corresponding  electric  currents  through  the  equivalent  internal 
resistance  Rint  and  the  leakage  resistance  Rieak,  respectively. 

With  considering  all  the  irreversibilities  discussed  in  the  above 
analysis,  the  power  output  of  the  fuel  cell  can  be  deduced  from  Eqs. 
(4)  and  (5)  as  [13] 


Pcell  =  -  AG  -  T0Stot  =  ^  (m  -  2^_m2)  (6) 

where 

d\  =  2  ne  sinh-1  (  )  +  2  ne  sinh-1  (  )-lnfl-J-) 

y2lo ,a  J  y2lo ,c  J  y  h,a  J 


1  engine 


Fig.  3.  The  schematic  diagram  of  an  endoreversible  Carnot  heat  engine. 


fuel  cell  can  be  defined  as 


hcell  '■ 


1  cell 


1  cell 


1  cell 


-AH  -( iA/neF)Ah°  -Ah 


y^-m-A2) 


(7) 


2.2.  An  endoreversible  Carnot  heat  engine 

For  a  practical  heat  engine,  there  are  invariably  thermal  resis¬ 
tances  between  the  working  substance  and  the  external  heat 
reservoirs.  For  the  sake  of  convenience,  the  heat  engine  in  the 
hybrid  system  is  assumed  to  be  an  endoreversible  Carnot  cycle  [38], 
as  shown  in  Fig.  3,  where  the  heat  Qh  is  transferred  from  the  high- 
temperature  heat  reservoir  (fuel  cell)  at  T to  the  working  fluid  at  T\ , 
and  Qi  from  the  working  fluid  at  T2  to  the  environment  at  T0.  When 
the  heat  transfer  between  the  reservoirs  and  the  working  fluid  is 
assumed  to  obey  Newton’s  law  [39],  we  have 

Qh  =  U1A1(T-T1)  (8) 

and 

Ql  —  U2A2(T2  —  T0)  (9) 

where  and  U2  denote  the  overall  heat-transfer  coefficients 
between  the  working  fluid  and  the  heat  reservoirs,  and  A\  and  A2 
are  the  heat-transfer  areas  between  the  working  fluid  and  the  heat 
reservoirs. 

As  illustrated  in  Fig.  2,  a  part  of  the  waste  heat  produced  in  the 
fuel  cell  is  directly  released  as  heat-leak  (Q/oss)  to  the  environment 
[7,39-42],  which  may  be  expressed  as  [43] 

Qtoss  =  KAi{T-To)  (10) 

where  I<  is  the  convective  and/or  conductive  heat-leak  coefficient, 
and  A/  denotes  the  effective  heat-transfer  area.  Combining  the  pre¬ 
vious  analysis  and  above  equations  yields 


■ln(1_G)  +  ^exp(&)’  = 

and  m  =  -Ah°  +  T  As°  +  RTln  (hAAl  )  -  RTdj 

(  Ph2o  J 

For  any  energy-conversion  device  that  takes  in  energy  to  pro¬ 
duce  work,  the  basic  definition  of  thermal  efficiency  is  the  useful 
work  output  divided  by  the  total  energy  input  [35-37].  Thus,  from 
the  thermodynamic  point  of  view,  the  efficiency  of  the  irreversible 


Qh  =  -AH  -  PceU  -  Qioss  =  —  AH  -  PceU  -  I<At(T  -  T0)  (11) 

According  to  the  definition  of  an  endoreversible  Carnot  cycle, 
one  has  Qh/Ti  =  Ql/T2  [39,44].  Using  the  above  equations,  it  can 
be  proved  that  when  A\ /A2  =  \/U2/U],  the  optimum  efficiency 
and  power  output  of  the  heat  engine  at  given  Qh  and  Ah  may  be 
expressed  as 

h  engine  —  1  —  Tq/(T  —  Qh/hTi) 


[T/T0  -  im2(l  -  rjceii )  +  m3(T/T0  -  1 )] 
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and 

P engine  =  Oh^I engine  =  <Tll  [if^ToU  -  Vcell)  ~  m3 (J  ~  To)] 


l  [T/T0-im2(l  -riceu)  +  m3(T/T0  -  1)]J 

where 

.  ,,  ,,  ./  /7T-  /77-\2  AAh°  I(A, 

m=AhUiU2/{VUi  +  VU2)  •  =  17,3  =  ^7 

and  =/\i  +A2  denote  the  overall  heat- transfer  area  of  the  heat 
engine  [38]. 

2.3.  An  ideal  regenerative  heat  exchanger 

As  illustrated  in  Fig.  1,  the  regenerator  in  the  hybrid  system 
works  as  a  heat  exchanger,  heating  the  inlet  reactants  from  the 
ambient  temperature  to  the  cell  temperature  by  using  the  high- 
temperature  outlet  gas  of  the  fuel  cell.  For  the  sake  of  simplicity, 
the  regeneration  process  is  assumed  to  be  ideal.  This  assumption 
is  reasonable,  because  the  efficiency  of  regenerators  with  the  val¬ 
ues  of  98-99%  have  already  been  reported  [45-47].  With  the  help 
of  perfect  regeneration,  the  fuel  cell  and  hence  the  whole  hybrid 
system  can  be  ensured  to  work  normally  and  continually  under  the 
condition  of  steady-state. 


2.4.  The  efficiency  and  power  output  of  the  hybrid  system 


Combining  Eqs.  (7)  and  (11 )— (13)  yields  the  following  expres¬ 
sions  of  the  efficiency  and  power  output  for  the  hybrid  system 
as 


h  hybrid  = 


P hybrid 

Qin 


Pcell  +  P( 


engine 


-AH 


—  hcell  + 


Pengine 

-AH 


and 


—  h  cell  + 


Qh  Pengine 
-AH  Qh 


=  riceu  +  [1  -  rjceu  -  {T/Tq  -  l)m3/(m2i)] 


x 


1 

[T/r0  -  im2(l  -  pceii)  +  m3{T/T0  -  1 )] 


} 


(14) 


P hybrid  ~  Qin  P  hybrid  —  Pcell  +  Pengine 
iA 


-(  neFAh°)  ('?ce,,+  I1  r?ceH  (t0  0 

1 


{ 


1  - 


m3  1 


[T/T0  -  im2(l  -  riceu)  +  m3(T/T0  - 


(m2i)J 

T»})  (151 


From  Eqs.  (14)  and  (15),  it  is  clearly  seen  that  the  efficiency  and 
power  output  of  the  hybrid  system  are  closely  dependent  on  the 
irreversible  losses  including  the  irreversibilities  within  the  fuel 
cell  itself  and  originating  from  the  heat  transfer  due  to  convec¬ 
tion/conduction  in  the  fuel  cell-heat  engine  hybrid  system.  In  the 
next  section,  numerical  predictions  will  be  studied  to  outline  how 
the  irreversible  model  based  on  the  above  analysis  can  provide  a 
valuable  tool  for  improving  the  system  performance. 


3.  General  performance  characteristics  and  optimal  criteria 

Optimal  design  and  analysis  of  a  system  require  a  thorough 
understanding  of  its  performance  limitations.  It  can  be  seen  from 
Eqs.  (14)  and  (15)  that  the  performance  of  the  fuel  cell-heat  engine 
hybrid  system  depends  on  a  set  of  thermodynamic  and  electro¬ 
chemical  parameters  such  as  the  operating  temperature  (T)  and  the 


Table  1 

Operating  conditions  and  performance-related  parameters  [20,48-52]. 


Parameter 

Value 

Operating  pressure,  p0  (atm) 

1 

Fuel  composition,  pHl ;  pH2o 

0.97;  0.03 

Air  composition,  po2 ;  Pn2 

0.21;  0.79 

Charge-transfer  coefficient,  ft 

0.5 

Number  of  electrons,  ne 

2 

Pre-factor  for  anode  exchange  current 

5.5  x  108 

density,  ya  (Am-2) 

Activation  energy  of  anode,  Eact,a 

1.0  x  105 

(J  mol-1 ) 

Pre-factor  for  cathode  exchange 

7.0  x  108 

current  density,  yc  (A  m  2 ) 

Activation  energy  of  cathode,  Eact,c 

1.2  x  105 

(J  mol-1 ) 

Electrolyte  thickness,  Le;  (p,m) 

20 

Activation  energy  of  O2-,  Ee[  (J  mol-1 ) 

8.0  x  104 

Pre-factor  of  O2-,  cr0  (S  m-1 ) 

3.6  x  107 

Ratio  of  the  internal  resistance  to  the 

1/100 

leakage  resistance,  k 

Anode  limiting  current  density, 

2.99  x  104 

(Am-2) 

Cathode  limiting  current  density,  ipc 

2.16  xlO4 

(Am-2) 

Faraday  constant,  F  (C  mol-1 ) 

96,485 

Universal  gas  constant,  R  (J  mol-1 1<-1 ) 

8.314 

Standard  molar  enthalpy  change  at 

-24,8303 

1073  K,  A h°  (J  mol-1 ) 

Standard  molar  entropy  change  at 

-55.5666 

1073  K,  As°  (J  mol-1  K-1 ) 

current  density  (2)  of  the  fuel  cell,  the  parameters  related  to  the  heat 
transfer  between  the  fuel  cell  and  the  heat  engine  and  the  heat-leak 
to  the  surroundings,  i.e.,  mi ,  m2 ,  and  m3.  Numerical  calculations 
are  performed  based  on  the  parameters  summarized  in  Table  1, 
which  are  derived  from  data  available  in  literatures  [20,48-52]. 
Their  values  are  kept  constant  unless  mentioned  specifically.  The 
fuel  composition  is  taken  as  97%  H2  +3%  FI20,  and  the  typical  oxy¬ 
gen  composition  in  the  ambient  air,  i.e.,  21%  02  +  79%  N2,  is  used  as 
oxidant.  Moreover,  the  enthalpy  and  entropy  of  reaction  are  gener¬ 
ally  not  strong  functions  of  temperature,  which  has  been  justified 
by  the  numerical  examples  [35,37]  for  hydrogen-oxygen  reaction. 
Thus,  we  can  invoke  the  assumption  that  the  changes  in  enthalpy 
and  entropy  across  the  reaction  are  independent  of  temperature. 

By  using  Eqs.  (14)  and  (15),  the  power  density  and  efficiency 
of  the  fuel  cell-heat  engine  hybrid  system  operating  at  a  tem¬ 
perature  range  of  1000-1500  K  and  a  current  density  range  of 
0-21600  A  m-2  are  presented  in  Figs.  4  and  5,  where  the  parameters 
m2  =  0.0001  m2  A-1  and  m3  =  0.001  are  chosen,  and  the  power  den¬ 
sity  Ptybrid  =  Phybrid/A.  It  can  be  clearly  seen  from  Figs.  4  and  5  that 
for  the  hybrid  system,  there  exist  a  maximum  power  output  Pmax 
and  a  maximum  efficiency  77max  as  the  current  density  i  and  oper¬ 
ating  temperature  T  are  varied.  Eqs.  (14)  and  (15)  show  clearly  that 
when  Ah0  is  assumed  to  be  independent  of  temperature,  there  is 
a  common  extremal  condition  of  dp  hybrid  I  dT  =  dPhybrid  I  dT  =  0  for  the 
hybrid  system.  It  implies  the  fact  that  there  is  a  common  optimum 
temperature  Topt  for  the  power  output  and  efficiency  of  the  system. 
This  characteristic  is  illustrated  more  clearly  in  Figs.  4b  and  5b.  Thus, 
in  the  practical  operation  of  the  fuel  cell-heat  engine  hybrid  sys¬ 
tem,  engineers  should  search  an  optimal  temperature  of  the  fuel  cell 
carefully  to  make  sure  that  the  whole  system  works  at  the  optimum 
condition. 

Eqs.  (14)  and  (15)  also  show  that  there  exist  two  extremal  con¬ 
ditions  of  dphybridldi  =  0  and  d Phybridldi  =  0  for  the  hybrid  system, 
but  dp  hybrid  I  di  =f=  dPhybrid/di .  It  means  that  the  current  density  at 
the  maximum  efficiency  {iv)  cannot  be  equal  to  that  at  the  maxi¬ 
mum  power  output  (ip),  as  shown  in  Fig.  6,  where  pm  and  Pm  are 
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Fig.  4.  The  power  output  density  as  a  function  of  the  current  density  and  temper¬ 
ature  of  the  fuel  cell,  where  the  parameters  m2  =0.0001  m2  A-1  and  m3  =0.001  are 
chosen.  The  values  of  other  parameters  are  listed  in  Table  1. 


the  efficiency  at  the  maximum  power  output  and  the  power  out¬ 
put  at  the  maximum  efficiency,  respectively.  It  is  seen  from  Fig.  6 
that  in  the  region  of  i<iv,  the  efficiency  and  power  output  of  the 
system  will  decrease  as  the  current  density  i  is  decreased;  while  in 
the  region  of  i  >  iP,  the  efficiency  and  power  output  of  the  system 
will  also  decrease  as  the  current  density  i  is  increased.  It  is  thus 
obvious  that  the  regions  of  i  <  iv  and  i  >  iP  are  not  optimal  from  the 
thermodynamic  point  of  view  although  the  hybrid  system  may  be 
operated  in  these  regions.  Therefore,  the  optimal  region  of  the  cur¬ 
rent  density  i  for  the  fuel  cell-heat  engine  hybrid  system  should 
be 


ir]  <  i  <  ip  (16) 

It  shows  that  iv  and  iP  are  two  important  parameters  of  the  hybrid 
system,  which  determine,  respectively,  the  upper  and  lower  bounds 
of  the  optimized  current  density.  In  the  practical  operation  of  the 
fuel  cell-heat  engine  hybrid  system,  engineers  should  choose  a  rea¬ 
sonable  current  density  according  to  Eq.  (16)  to  ensure  that  the 
system  is  operated  at  the  optimal  region. 

According  to  the  optimum  criterion  of  the  current  density  and 
Fig.  6,  one  can  further  determine  the  optimum  regions  for  the  effi¬ 
ciency  and  power  output  as 

Pm  <  Phybrid  —  ^max  (17) 
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Fig.  5.  The  efficiency  as  a  function  of  the  current  density  and  temperature  of  the 
fuel  cell.  The  values  of  the  relevant  parameters  are  the  same  as  those  used  in  Fig.  4. 


and 

l)m  <  1J hybrid  —  *7max  (18) 

When  the  hybrid  system  is  operated  in  these  optimum  regions, 
the  power  output  will  increase  as  the  efficiency  is  decreased,  and 
vice  versa.  The  above  results  show  that  Pm ax,  Vmax,  Pm,  and  rjm  are 
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Fig.  6.  The  curves  of  efficiency  and  power  density  varying  with  the  current  den¬ 
sity  of  the  fuel  cell-heat  engine  hybrid  system  for  the  parameters  T  =  1183  K  and 
m2  =  0.00063  m2  A-1,  where  iv  and  iP  are  the  current  densities  at  the  maximum 
efficiency  ^max  and  maximum  power  output  density  P^ax,  respectively.  The  values 
of  other  parameters  are  the  same  as  those  used  in  Fig.  4. 


Y.  Zhao,J.  Chen  /  Journal  of  Power  Sources  186  (2009)  96-103 


101 


also  four  important  parameters  of  the  fuel  cell-heat  engine  hybrid 
system.  Pmax  and  rjmax  determine  the  upper  bounds  of  the  power 
output  and  efficiency  of  the  system,  while  Pm  and  rjm  determine  the 
allowable  optimum  values  of  the  lower  bounds  of  the  power  output 
and  efficiency.  Obviously,  the  four  important  parameters  depend 
closely  on  some  of  the  system  parameters  and  can  be  numerically 
calculated  for  the  given  values  of  these  parameters. 

4.  Results  and  discussion 

It  is  noteworthy  that  the  developed  system  model  can  be  used 
to  investigate  the  effect  of  the  various  operating  and  design  param¬ 
eters  on  the  performance  of  the  hybrid  system.  In  this  section,  a 
parametric  analysis  is  carried  out  based  on  the  performance  cri¬ 
teria  for  the  variation  of  operating  conditions  to  characterize  the 
system  behavior. 

4  A.  Effect  of  and  m 3 

By  using  Eqs.  (14)  and  (15),  it  can  be  easily  proved  that  the 
power  output  and  efficiency  of  the  hybrid  system  are  monotoni- 
cally  increasing  functions  of  mi.  Its  physical  meaning  is  quite  clear. 
Because  mi  is  a  parameter  to  measure  the  irreversibility  of  finite- 
rate  heat  transfer  in  the  heat  engine,  the  larger  the  parameter  mi , 
the  smaller  the  heat  transfer  irreversibility  in  the  heat  engine  and 
the  better  the  performance  of  the  hybrid  system. 

From  Eqs.  (14)  and  (15),  it  can  be  also  proved  that  the  power 
output  and  efficiency  of  the  hybrid  system  are  monotonically 
decreasing  functions  of  m3.  Its  physical  meaning  may  be  easily 
explained  as  follows:  m3  is  a  synthesized  parameter  to  measure  the 
heat-leak  irreversibility  from  the  fuel  cell  to  the  surroundings  and 
the  heat  transfer  irreversibility  in  the  heat  engine.  The  smaller  the 
parameter  I<At  is  and  the  larger  the  parameter  mi  is,  i.e.,  the  smaller 
the  heat-leak  irreversibility  of  the  fuel  cell  and  the  heat  transfer 
irreversibility  in  the  heat  engine  are,  the  smaller  the  parameter 
m3  is,  and  consequently,  the  better  the  performance  of  the  hybrid 
system  will  be. 

4.2.  Effect  of  m2 

The  parameter  m2  is  a  colligation  measurement  for  the  systemic 
structure.  By  using  Eqs.  (14)  and  (15),  it  may  be  proved  that  there 
is  a  common  extremal  condition  of  drjhybridldrn2  =  dPhybridl^m 2  =  0 
for  the  hybrid  system.  It  implies  the  fact  that  for  given  values  of 
the  current  density  i  and  temperature  T  of  the  fuel  cell,  there  is 
a  common  optimum  (m2)0pt  at  which  both  the  power  output  and 
efficiency  attain  their  maxima,  as  shown  in  Fig.  7,  where  (m2)p  and 
(m2)p  are  the  values  of  m2  when  the  current  density  is  equal  to  in 
and  iPt  respectively.  According  to  Eq.  (16),  we  can  further  determine 
the  optimal  region  of  the  parameter  m2  as 

(m2)P  <  (m2)opt  <  (m2)r,  (19) 

Obviously,  (m2)r,  and  ( m2)p  are  two  important  parameters  of  the 
hybrid  system,  which  determine,  respectively,  the  upper  and  lower 
bounds  of  the  optimized  parameter  (m2)opt- 

4.3.  Effect  of  m3/m2 

When  m\  oc,  m2  =  0  and  m3  =  0,  but 
m3/m2  ={neFI<AiT0)l(-A  Ah°)  is  a  parameter  which  is  inde¬ 
pendent  of  m\.  In  fact,  m3/m2  is  only  a  parameter  to  measure  the 
performance  of  the  fuel  cell  and  is  independent  of  the  heat  engine. 
Figs.  8  and  9  show  clearly  the  effect  of  m3\m2  on  the  power  output 
and  efficiency  of  the  system  under  the  condition  of  m\  -*  oo.  It 


Fig.  7.  The  curves  of  the  efficiency  and  power  density  of  the  hybrid  system  varying 
with  m2  for  the  parameters  7=  1183  K,  m3  =0.001,  where  (7712)77  and  (m2)P  are  the 
values  of  m2  when  the  current  density  is  equal  to  and  iP,  respectively.  The  values 
of  other  parameters  are  the  same  as  those  used  in  Fig.  4. 


Fig.  8.  The  curves  of  (a)  the  efficiency  and  (b)  power  output  density  of  the  hybrid  sys¬ 
tem  varying  with  current  density  of  the  fuel  cell  for  differently  given  values  of  m3  /m2 , 
where  parameters  m2  =  0,  m3  =  0,  7=  1183  K,  and  the  values  of  other  parameters  are 
the  same  as  those  used  in  Fig.  4. 
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Fig.  9.  The  curves  of  the  efficiency  and  power  density  of  the  hybrid  system  varying 
with  the  fuel  cell  temperature  for  differently  given  values  of  m3/m2,  where  param¬ 
eters  m2  =  0,m3=  0,  i  =  2000  A  m~2 ,  and  the  values  of  other  parameters  are  the  same 
as  those  used  in  Fig.  4. 


is  observed  from  Figs.  8  and  9  that  both  the  power  output  and 
efficiency  of  the  system  will  increase  obviously  as  the  value  of 
m3/m2  is  decreased.  The  reason  is  quite  obvious.  The  smaller 
the  heat-leak  loss  from  the  fuel  cell  to  the  surroundings  and  the 
larger  the  enthalpy  change,  the  smaller  the  parameter  m3/m2,  and 
consequently,  the  better  the  performance  of  the  hybrid  system. 


4.4.  Several  interesting  cases 


(1)  When  the  influence  of  the  heat  transfer  irreversibility  in  the 
heat  engine  is  negligible,  i.e.,  the  overall  heat-transfer  coeffi¬ 
cients  and  U2  tend  to  infinity,  rrq  -*  oo,  m2  =  0,  m3  =  0,  and 
Eqs.  (12),  (14)  and  (15)  may  be,  respectively,  simplified  as 


i  To 

4 engine  ~  *  j  —  4c 


4  hybrid  ~  4 cell  + 

and 


4c 


^hybrid  =  (“ 


(->•){ 


4  cell  + 


1  4  cell 


(?)(£ 


(20) 

(21) 


where  r]C  is  the  Carnot  efficiency  and  m4  =  m3/m2  = 
( neFKAiT0)l(-AAh° ).  In  such  a  case,  there  still  exist  the 
maximum  values  for  the  power  output  and  efficiency  of  the 
hybrid  system  and  the  optimal  values  of  the  current  density  and 
operating  temperature  of  the  fuel  cell,  as  shown  in  Figs.  8  and  9. 
Flowever,  in  a  practical  fuel  cell-heat  engine  hybrid  system, 
the  overall  heat-transfer  coefficients  and  U2  are  always 
finite.  m\  cannot  tend  to  infinity  and  m3  is  always  larger  than 
zero.  Thus,  the  results  obtained  by  using  Eqs.  (14)  and  (15) 
to  analyze  the  performance  of  a  fuel  cell-heat  engine  hybrid 
system  will  have  more  realistic  meaning  for  the  optimal  design 
of  the  system. 

(2)  When  the  heat-leak  from  the  fuel  cell  to  the  environment  is 
negligible,  i.e.,  KA{  0,  m3  =  0,  and  Eqs.  (12),  (14)  and  (15)  may 
be  simplified  as 


4  engine  —  1 

1 

(23) 

[T/To  -  im2(l  -  r]cen)] 

4  hybrid  =  1 

1  —  ^Icell 

(24) 

T/Tq  -  im2(l  -  4 cell) 

and 


^hybrid  —  (  nep  Ah°) 


_ 1  4  cell _ 

T/Tq  -  im2(l  -  rjceu)_ 


(25) 


(3 )  When  both  the  irreversibilities  due  to  the  heat  transfer  between 
the  fuel  cell  and  the  heat  engine  and  the  heat-leak  from  the 
fuel  cell  to  the  environment  are  negligible,  i.e.,  m2  =  0,  m3  =  0, 
and  m3\m2  =  0,  we  can  obtain  a  simple  system  model,  and  con¬ 
sequently,  the  results  obtained  above  can  be  simplified.  For 
example,  Eqs.  (14)  and  (15)  can  be,  respectively,  rewritten  as 


4hybrid  —  4cell  +  0  4cell)4c  (26) 

and 


P hybrid  —  ^  ^  p  ^  [*7 ce// +  0  4cell)4c]  (27) 

As  shown  clearly  by  the  solid  curves  in  Figs.  8  and  9,  the  heat- 
leak  does  indeed  have  an  obvious  effect  on  the  power  output 
and  efficiency  of  the  hybrid  system. 

(4)  When  the  fuel  cell  of  the  hybrid  system  is  totally  irreversible, 
it  means  that  the  fuel  cell  works  at  the  open  (or  short)  circuit 
condition.  In  such  a  case,  Pcen  =  0  and  i)cen  =  0,  but  the  heat  engine 
still  works  normally,  and  the  hybrid  system  has  an  efficiency 
and  power  output  represented  by  the  following  equations: 


[i  (T  i 

\  ™3  I 

[  U 

)  ( m2i ). 

1  - 
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[T/To  -im2  +  m3(T/T0  - 


T»} 


(28) 


and 


r hybrid  —  ^engine  : 


1  - 

(1-1 

\  in  3  - 

\  neF  ) 

\T0 

)  ( m2i ). 

x  n- 


[T/T0  -  im2  +  m3(T/T0  - 


to) 


(29) 


(5)  When  the  heat  engine  of  the  hybrid  system  is  totally  irre¬ 
versible,  P engine =  and  consequently,  the  hybrid  system 

reduces  to  the  simple  system  of  the  pure  fuel  cell  and  its  optimal 
performance  has  been  given  in  Ref.  [13]. 


To  sum  up,  by  comparing  the  results  obtained  by  the  hybrid 
system  with  those  derived  from  the  simple  systems  only  with  the 
fuel  cell  or  the  heat  engine,  it  is  easily  seen  that  the  performance  of 
the  hybrid  system  is  much  better  than  that  of  the  simple  systems. 
It  indicates  that  the  fuel  cell-heat  engine  hybrid  power  system  is 
an  obvious  choice  for  the  practical  applications. 


5.  Conclusions 


The  present  paper  shows  a  concept  of  integrating  a  heat  engine 
into  an  existing  fuel  cell  model.  It  focuses  on  the  system  modeling 
of  the  fuel  cell-heat  engine  hybrid  and  tries  to  improve  the  whole 
system  performance  through  numerical  simulations.  The  various 
irreversible  thermodynamic  and  electrochemical  losses  are  iden¬ 
tified,  and  the  individual  contribution  of  these  irreversibilites  to 
the  whole  system  performance  is  investigated  through  parametric 
study.  The  optimum  criteria  of  some  important  parameters  such 
as  the  power  output  and  efficiency  of  the  system,  the  current  den¬ 
sity  and  operating  temperature  of  the  fuel  cell,  and  the  parameters 
related  to  the  systemic  structure  are  given.  Consequently,  the  per¬ 
formance  characteristics  of  the  hybrid  system  are  described  and 
optimized  from  a  thermodynamic  point  of  view.  This  new  method 
may  be  easily  extended  to  other  fuel  cell  hybrid  systems  to  develop 
irreversible  models  suitable  for  the  optimal  energy-management 
strategies  of  fuel  cell  hybrids. 
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